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We consider N identical two-level atoms coupled to an optical cavity, which is coherently driven by 
an external field. In the limit of small atomic excitation, the reflection and transmission coefficients 
for both fields and intensities are calculated analytically. In addition, the frequency content of 
the cavity field and hence also the emission spectrum of the cavity is determined. It is discussed in 
particular how individual collisional dephasing and common atomic energy-level fluctuations prevent 
the cavity field from being in a coherent state, which in turn affects the outgoing fields. 
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I. INTRODUCTION 



In the field of cavity-quantum-electrodynamics, the 
coupling between a single two-level system and the 
electro-magnetic field is enhanced by a cavity, which al- 
lows for the study of light-matter interactions in the most 
fundamental waji Such a setup has been implemented 
in various physical systems, e.g.: A Rydberg atom cou- 
pled to a micro-wave cavity^, an alkali atom coupled to 
an optical cavity^, a super-conducting qubit coupled to 
a transmission line resonator— and a semi-conductor 
quantum dot coupled to a photonic-crystal cavity^. En- 
sembles of two-level atoms, each coupled weakly to the 
cavity field, are also able to present an effective strong co- 
operative couplingiSr— . The light-matter coupling gives 
rise to a normal-mode splitting of the resonance frequen- 
cies [known as the vacuum-Rabi splitting as explained 
by the Jaynes-Cummings models for a single atom], 
which can be detected dynamically as oscillations be- 
tween atomic and optical excitations^, as a double-peak 
in the steady-state fluorescence spectrum^, or as a dou- 
blet structure in the cavity-transmission profile- 5 . The 
transmitted and reflected fields can be exploited in real- 
time detection and control of trapped atom s 16 ' 17 , and the 
reflection and transmission profile of these fields has pre- 
viously been calculated using both classical^ 5 - and quan- 
tum theories^. While the atomic-dipole-moment decay 
was most likely lifetime limited in the traditional ex- 
periments with free atoms in vacuum^— , the modern 
solid-state implementation has called for an increased 
attention to dephasing mechanisms in the case of both 
fluorescence^ and transmission measurements^. The 
present work derives analytical expressions for the trans- 
mitted and reflected intensities, which under the influ- 
ence of atomic dephasing noise do not coincide with the 
modulus square of the transmitted and reflected fields. 
In addition, the spectrum of the field emitted from the 
cavity is determined. 
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FIG. 1. The schematics of the physical system under consid- 
eration. An optical cavity contains a cavity field, a c , which 
couples equally to N two-level atoms. The field-decay rates 
through the two mirrors are denoted by Ki and K2, which 
connect the cavity field to the coherent input field, /3, the 
reflected field, or, and the transmitted field, &t- 

II. AN EMPTY CAVITY WITH PHASE NOISE 

Before considering the real problem of a cavity con- 
taining N atoms subjected to various decay processes, 
we wish to motivate the efforts by a simpler example: 
An empty cavity subjected to phase noise and driven co- 
herently by the field /3 in = /3e _I " L * (Fig. Q] without atoms, 
N = 0). In the frame rotating at the driving frequency, 
wl, the Hamiltonian of this simple system is given by: 

H = hA c ala c + ihy/^iP&l - /3*a\), (1) 

where A c = uj c — ljl is the detuning from the cavity res- 
onance frequency, uj c . The input field amplitude, /3, is 
normalized such that |/3| 2 is the number of photons inci- 
dent on the cavity input mirror per second, and the field- 
decay rate, k 1; at the input mirror connects the cavity 
field, a c , to the incident and reflected fields, /? and ajj,, by 
the input-output formalism of Collett and Gardiner—: 

a R = \/2K^a, c - (fl + vi), 

nz — „ „ ( 2 ) 
ax = v ZK2a c — V2- 

The second line determines the transmitted field, St, us- 
ing the field-decay rate, K2, of the output mirror. In the 
above equations v\ and 62 operate on vacuum states and 
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will have no further impact on the calculations. In order 
to account for decay processes, we take the approach of 
the master equation, || = jf-[H, p] + C(p), where the 
Lindblad part, C(p), is a sum of terms: 



Cm(p) 



'T.^p'rnPraP + pCm^m) + CmpC^ 



(3) 



The leakage of photons from the cavity is modeled by a 
Lindblad term with C\ = \f2ka ci where k — Ki + k 2 is 
the total cavity-field decay rate. 

Now, we wish to include phase noise in the cavity. 
Physically, this could arise from a fast jitter of one of 
the mirrors such that the cavity resonance frequency be- 
comes u' c = lj c + s(t), where e(t) is a real function ful- 
filling (e(t)} = 0, and (s(t)e(t')) = ^5(t - t'). Com- 
putationally, the time-dependent term, H^t = fe(i)aja c , 
could be added to the Hamiltonian; however, the fact 
that the noise is delta-correlated enables a simpler, time- 
independent modeling using a Lindblad term with C*2 — 
y^^-dja c , see App.lA"! The master equation now leads to 

the dynamical equations for the cavity field and photon 
number: 



9< 



dt 
djajac) 
dt 

where F = 



= -0^A c )<a c )+x/2^/3, (4) 
= -2K(ala c ) + ^hT 1 (p(al)+P*{a c )), (5) 



l 



In steady state, the mean field am- 



plitude and photon number can be written: 



T + iAr 



(ala c ) 



(6) 
(7) 



This result is important for calculating the reflected S^Sr 
or transmitted d T dx intensity using Eqs. @, and it is 
evident that the introduction of phase noise (rjn < oo) 
prevents the cavity-field steady state from being coher- 
ent, which would otherwise lead to the simpler relation, 
(dja c ) = | (d c ) | 2 . The gradual change from a coher- 
ent state, when the contribution (KTjit) -1 to Eq. ([7]) in- 
creases from zero, is shown by the Wigner function^ in 
Fig. [5] In our real iV-atom problem to be discussed in the 
next section, the atoms will effectively add phase noise 
to the cavity field, and, similar to the present example, 
the transmitted and reflected intensities cannot simply be 
calculated as the modulus square of the field amplitudes. 
We note that the impact of phase noise in a cavity on the 
density matrix itself has been studied experimentally by 
Wang et al.—. 

III. A CAVITY FILLED WITH ATOMS 

Now, consider the entire setup of Fig. [T]with N identi- 
cal two-level atoms placed in the coherently driven opti- 
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FIG. 2. The Wigner function for the cavity field for various 
contributions of phase noise. In all panels, (a c ) = 2, and the 
value of (fCTjit) -1 is varied as: (a) 0.0, (b) 0.1, (c) 0.3, and (d) 
1.0. 



cal cavity. In addition to the empty-cavity Hamiltonian 
of Eq. ([T]), we add the contribution from atoms: 



TiA., 



N 



JV 

Tig^f^ + a^al). (8) 



Here A a = w a - wl is the detuning of the driving fre- 
quency, wl, from the atomic resonance frequency, w a . 
The coupling constant, g, is taken to be equal for all 
atoms, the j'th of which is described by the Pauli op- 



We restrict the cal- 



(i) 

erators, with k 
culations to the regime of low excitation probability, 

JW = \{{v^J +1)<1- This allows for the Holstein- 
Primakoff approximation, in which each atom is modeled 
by a harmonic oscillator: 



a 



GO 



at, a ( z j) 2a] aj - 1. 



(9) 



With this approximation, the entire system Hamiltonian 
reads: 



H 



hA c a' c a c 



N 



N 

KA a 22 



+ Tig VVdJa c + djdj) + ih^/2K^(f3al - f3*al), 



3=1 



(10) 



GO 



where the constant-energy shift from the — 1-term in a 
in © has been omitted. We assume that the phase noise 
from the jittering mirror (Sec. [TTJ) is absent (rjit = oo, 
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i.e. C2 = 0), but the cavity leakage is retained and vari- 
ous atomic decay processes will be included. Firstly, for 
each two-level atom the decay of population with rate 7y 



preted in the Schrodinger picture by the Monte-Carlo- 
wave- function approach 23 . This method states that the 
wave function should evolve under the non-Hermitian 



is modeled by C3 = ^/7jjcr_ 



where the ar- Hamiltonian, H - ^- 7j < 



row denotes the replacement suggested by Eq. (J9J> . Sec- 
ondly, dephasing of each individual atomic dipole mo- 
ment is accounted for by C4 = ~~ * \f^^j®j> 
where r is the dephasing rate^ 3 . [the constant term —1 
in ^ has been removed as it makes no contribution to 
the Lindblad operator Q when C m is hermitian] . Note 
that this Lindblad term can be derived as the average 
effect of fast random fluctuations of the resonance fre- 
quency of each atom along the lines of App. Due to 
the Holstein-Primakoff approximation, the above Lind- 
blad terms, C3 and C4, for atomic-population decay and 
dipole-dephasing correspond exactly to the energy and 
phase decay in the harmonic oscillator—. Finally, in 
addition to individual atomic decay and dephasing pro- 
cesses, we add the possibility of dephasing mechanisms 
common to all atoms. For instance, a fluctuating mag- 
netic field could vary the energy levels of an ensemble of 
magnetic dipoles in a uniform manner, and laser-intensity 
fluctuations could cause common energy-level variations 
of atoms in an optical dipole trap. This scenario is mod- 
eled by an additional fluctuating term in the Hamilto- 
nian, Hi — he(t) ! where we assume the correla- 
tion time of fluctuations to be fast compared to any other 
time scale: (e(t)) = 0, and (e(t)e(t')) = ^S(t - t'). As 
explained in App. such a fluctuating term can be cast 



into Lindblad form ([3]) with C5 = 1/ 4 E 



Based on the results in Sec. [TT1 we expect atomic 
phase noise to destroy the coherence of the cavity field. 
Conversely, in absence of phase noise (r = r' = 00) 
but maintaining population decay (711 > 0), the cavity 



field is indeed coherent with (aja c ) 



provided 



that the Holstein-Primakoff approximation holds. To see 
this, observe that the Hamiltonian (fTUf is on the form 
H = k qj^o^ak, where j, k are indices for atoms, the 
cavity, and the external (c-number) driving field, and qj^ 
are time-independent complex numbers. In the Heisen- 
berg picture: gjdj(i) = jf^J2k1j,kak(t), which in turn 
means that the vector a of all annihilation operators, a.j, 
evolve as a(t) = exp(— iQt)a(0), where Q is the ma- 
trix of entries q^k- Now, if at time t = all oscillators 
are in a coherent state with amplitudes a (e.g. all atoms 
and the cavity in the vacuum state and the external 
driving field in a coherent state with amplitude 0) such 
that a(0) — a\ip), we find at later times: a(t) — 
exp(— j-Qt)a i.e. the state remains an eigenstate for 
all annihilation operators, and in particular the steady 
state of each oscillator must be coherent. Of course, 
this observation holds in the Schrodinger picture as well. 
Now, if population decay (with rate 7^) of an oscillator 
is introduced, i.e. the Lindblad term with C — *J^fjO,j is 
included in the calculations, the evolution can be inter- 



and occasionally be sub- 
jected to a quantum jump, \ip) — > cij \tp) = otj \ip}- The 
latter clearly preserves \ifj), and the non-Hermitian term 
is of the form discussed above and thus also preserves the 
coherent state. 



IV. STEADY-STATE SOLUTIONS 



A. Exact expressions 



In order to calculate the relevant physical variables in 
steady state, the dynamical equations of the cavity field 
and atomic dipoles are deduced as a first step: 



d(a c 
dt 

d {dj 



N 



= -(k + iA c ) (fie) - igJ2 <%> + V^iX ( u ) 



-(71 +«A a ) (a,j) - ig(a c ) , 



(12) 



where 7^ = 7 + 77 + • The steady state for the atomic 



operators becomes (cij) = , 
the cavity-field steady-state value 



which in turn leads to 



g 2 N 



(k-HA c )(1 + u)' (k-MA c )( 7 _l -HA a ) 

(13) 

The complex field-reflection and transmission coefficients 
can be immediately deduced using Eq. ([2]): 



(aR> 



(or) 



2ki 



(k + iA c )(l + v) 
(k + iA c )(l + v)' 



- 1 



(14) 
(15) 



These equations incorporate the dephasing mechanism 
through the definition of 7^ and have been presented 
previousl y 15 ' 18 . However, a correct treatment of intensi- 
ties requires the computation of the number of photons 
in the cavity. To this end, we deduce the following equa- 
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tions of motions: 



2k (ata c ) + [P (oj) + /?* (a c )] 



N 

l 9^2 - (oca] 



(16) 



3=1 



— (aja,-) = - (k + 7j_ + iA ac ) (aja,-) + \f2k~lf3* (aj) 



N 



ig(ala c ) +ig^2 



fc=i 



9 / At . 
— ( a o,- 
dt \ k 3 



= + ig 



(&%) - (ac4)] 



2 ( 4% 



4,i " 1 



21 

2 



(17) 



(18) 



where A ac = A a — A c has been defined. These equations 
can be solved in steady state [the details are presented 
in App. |B]: 



(&la c ) = 



2t' 



(19) 



2g 2 l(a c )| 2 ^ 

7l + A 2 



JV-yii 



7X--D 



21 

2 



where Z? is a constant given by: 
D=2 Kll] [(K + 1± ) 2 + Al c ] 



+ Ang 2 {n + 11 _) 



1 + f 7||r 
1+37117- 



2k 



(20) 



(21) 



Equations (TT9")) and (f2"0)) present one of the main re- 
sults of the present paper; they determine the relevant 
properties of the coupled atom-cavity system, and with 
these at hand the intensity-reflection coefficient, R = 

(IrSr) /\f3\ 2 , and the intensity-transmission coefficient, 



T 



a T ) /|/3| 2 , can be stated using Eq. ©: 



R=\r\ 



T=\t\ 



Ak\ 



k 2 + A 2 c |l + tf 

{a\a c ) 




(22) 



(23) 



Here r and t are the reflection and transmission coef- 
ficients, respectively, for the fields given by Eqs. ([14]) 
and (fT5|) . and the ratio (aja c ) /| (a c ) | 2 is given by 
Eq. ([T9"l) . In the absence of phase noise, this ratio is 
unity and the above expressions simplify into R = \r\ 2 
and T — \t\ 2 . In order to exemplify the practical impact 
of phase noise, consider Fig. [3l which shows the trans- 
mission profile for various decay parameters. In com- 
parison to the case with no phase noise (dotted line) , the 




FIG. 3. The cavity intensity-transmission coefficient, T, as 
a function of driving frequency, ui-l. For all curves, g = 2, 
N — 5, 2ki = 2^2 = K = 1, w a = u) c = 0, and 711 = 2. Dotted 
line: No phase noise (70. = dashed line: Only common 
phase noise (7_i_ = 2 7 ||), solid line: Only individual phase 
noise (7_i_ = 27||). The open circles show the modulus square 
of the field-transmission coefficient, \t\ 2 , when 7j_ = 2711. 



transmission profile broadens and decreases in magnitude 
when phase noise is added (solid and dashed lines) since 
j± increases from to 27y, which in turn decreases the 
magnitude of \t\ 2 in Eq. (f23|) . However, the transmission 
remains higher than \t\ 2 (open circles) due to the effect of 
Eq. (|19p . and this effect is larger for the common phase 
noise (dashed line) as compared to individual case (solid 
line). On a relative scale, this effect is most pronounced 
at ojl = [at atomic resonance, A a = 0]. 



B. Approximate results, limits 



Equations (|19p and (|2"0"|) are quite involved, but they re- 
duce to simpler forms in interesting limiting cases. Con- 
sider the fraction (which is common in both equations): 



7± 



JV-yjl 
2t' 




when t = t' = 
~ l when N = 1, 
when t = 00, 
when t = 00. 



00, 



(24) 



This shows, that in absence of phase noise (r = t 1 = 00), 
the square brackets in Eqs. (fT9"|) and (f2"01) reduce to 
unity. This confirms the coherent state of the cavity 
field, (aja c ) = | (a c ) | 2 , and Eq. ([20]) reduces to the semi- 
classical expression for the atomic excitation probability 
in a classical field [with Rabi frequency, x — %g (a c )]- In 
the case of one atom, N = 1, there is no difference be- 
tween common and individual phase noise, and the two 
decay channels can be considered effectively as a single 
channel with characteristic rate, r _1 + (t') _1 . If only 
one or the other type of dephasing is present [individ- 
ual collision-like dephasing only (r' = 00, r < co) or 
common dephasing only (r = 00, r' < 00)], its effect is 
proportional with the rate of the process, and for com- 
mon phase noise it is enhanced by a factor N. If we 
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FIG. 4. Both panels show on the vertical scale the 
height, h, of the Lorentzian contribution: (aja c ) = [1 + 

2 

h- 7 2~!^a 2 ] ' I I 2 ' wrme the atomic decay parameters are 
varied. Solid line: Only individual phase noise (r' = oo). 
Dashed line: Only common phase noise (r = oo). The dotted 
line shows the cooperativity parameter C = f^r~~- In panel 
(a) the dephasing time (r or r') is varied with y\\ = 2. In 
panel (b) 711 is varied with the dephasing time (r or r') set 
to unity. For both panels: g = 2, N = 5, 2ki = 2k2 = K — 1, 
and A ac = 0. 



7^ 



tI+a! j 



the 



write Eq. (JTSJ) as (aja c ) = [I + /i • 
pre-factor /i of the Lorentzian function in Eq. (fTI))) fulfills: 



ft < 



2C 

■ 711 



c 



always, 
when i, 



. < 21 
^ 2 • 



(25) 



where C = i s the cooperativity parameter. The 
maximum value C is obtained in the asymptotic limit 
7|| — > 0. The variations in h as a function of the atomic 
decay parameters is illuminated further in Fig. 2) In 
both panels, the cooperativity parameter C (dotted line) 



clearly presents an upper limit for h. Also, for these 
examples the common phase noise (dashed line) gives a 
larger contribution to h than individual phase noise (solid 
line) as indicated by Eq. (|25p , but otherwise present qual- 
itatively similar features. In fact, when the decay is life- 
time dominated (-j- - , -p ; right-hand-side limit of the 
two graphs) the height h decreases as (dephasing time) -1 
in panel (a) and as j^ 2 in panel (b) as predicted by 
Eq. (|23|) [note, C scales as y^ 1 in this limit]. Conversely, 
if the dephasing mechanisms dominate the lifetime de- 
cay (-J- <C 7, 77; left-hand-side limit of the two graphs), 
the cooperativity parameter C is the important figure of 
merit. In the limit 711 — > 0, one obtains h — > C as ex- 
emplified in Fig. HJb). This limit is reached whenever 



<C mm 



kN- 1 , (Nt)- 



( k+Xl )(1+- 



A2 



For fi- 



nite 7|| = 2, as exemplified in Fig. Eta), an approximate 
linear scaling of h with C is seen to apply. 



V. THE CAVITY EMISSION SPECTRUM 

The cavity field, a c , may possess frequency components 
different from the driving frequency, cjl, when phase 
noise is present. The spectral density characterizes this 
effect: 

1 f°° 

Sa c (w) = 5- / (at(t)a c (t + r)) e lA ^r, (26) 

where A = uj — cjl accounts for the rotating-frame pic- 
ture of a c . Physically, Ss. c (w)du) measures how much op- 
tical energy is present in the cavity within a frequency 
bandwidth, Sui, around u>. We note that J^ o S'a ( .(w)(iw = 
(qt& c ). With the computational details presented in 
App. \C\ the spectrum reads for an empty cavity sub- 
jected to phase noise (as described in Sec. HI]) : 



t/tt 



T 2 + (UJ C - u) 2 KTjit 



+ 8(uj - cjl) 



( 27 ) 

where | (a c ) | 2 refers to the steady-state cavity field de- 
pending on cjl through Eq. ([6]). The spectrum is divided 
into two parts: A broad-band, incoherent Lorentzian 
term which is only present when the cavity is subjected 
to phase noise (rjit < 00), and a coherent part oscillat- 
ing exactly at the driving frequency, cjl- When atoms are 
included (as described in Sec. the spectrum becomes: 



(k + 7±)(«7_l +9 2 N) + 



ir\g 2 N +(k + i[u c - u])(-y± + i[cJ a - uj})\ 2 7^ + (w a - cj l ) : 



+ S(w — u>l) 



(28) 



where h is the peak height discussed around Eq. (1251) . again, the spectrum is divided into a broad-band, in- 
We remind that h = in absence of phase noise. Once 
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FIG. 5. 



The cavity emission spectrum [note the log scale] 
modeled with a resolution bandwidth of 7tk d = 1(T 2 . For all 
curves, g — 2, TV = 5, 2ki = 2^2 = K = 1, w a = w c = 0, and 
7|i = 2. Solid and dashed lines: Only common phase noise 
(7_l = 2711) with driving frequencies, ojl = and ojl = 8, 
respectively. Dotted line: No phase noise, ul = 0. 



coherent part and a delta-function term accounting for 
the coherent part. In both Eqs. ([27]) and (|28]l. the left- 
most fraction in the square brackets is a function with 
unity area being proportional to T(uj) — |i(oj)| 2 , where 
T and t are the intensity and field-transmission profiles, 
but with the driving frequency, ojl, replaced by the ob- 
servation frequency, uj. These fractions are then multi- 



plied by the steady-state value of (aja c ) 



(the 



right-most fractions in the two equations together with 
I (a c ) | 2 ), which thus becomes the area of the incoherent 
part of the spectrum. Due to the delta-functions, the co- 
herent part contributes the remaining area, | (d c ) | 2 , such 
that the total area amounts to (djd c ) as expected. Ex- 
amples of spectra are presented in Fig. [5] with parameters 
closely resembling the dashed-line intensity-transmission 
profile of Fig. [3] As can be seen, it is only the magnitude 
of the incoherent spectrum which depends on the driv- 
ing frequency, wl; its shape corresponds to the difference 
between the dashed line and the circles of Fig. [3] The 
sharp peaks resemble the coherent part, which follows 
wl- While the coherent part of the reflected and trans- 
mitted fields is related to the cavity field by the input- 
output relations Eq. @, the incoherent, broad-band part 
of the spectrum is simply 2/ti or 2k,2 times the incoher- 
ent cavity-field spectrum, Sa (ui), for the reflected and 
transmitted fields, respectively. We stress that the cavity 
emission spectrum [through its exit mirrors] depends on 
the field-correlation function, (dj(t)d c (t + r)), and does 
not coincide with the atomic fluorescence spectrum^, 
which is radiated into all space and depending on the 
atomic-dipole correlation function, (a+ (t)&- (t + r)) . 



VI. DISCUSSION 

The previous sections have identified the practical in- 
fluence of atomic phase noise in two cases. Firstly, 



in a transmission or reflection measurement where the 
frequency, ojl, of the coherent driving field is varied, 
the intensity-transmission and reflection profiles, T(wl) 
and -R(wl), do not coincide with the modulus square 
of field-transmission and reflection profiles, |£(wl)| 2 and 
|r(wL)| 2 , according to Eqs. (|2"2"j) and (|2U)) . As a result, 
the pure dephasing-part of the atomic decay processes 
can be identified in experiment by e.g. comparing T and 
|i| 2 , since the ratio should present a Lorentzian feature, 

|7p = 1 + 7 a _|7^2 ; provided that the cooperativity pa- 
rameter is not negligible. The transmitted intensity can 
be deduced simply by using a photo-detector, whereas 
the value of \t\ 2 could be determined from measuring two 
orthogonal quadratures of the transmitted field in a ho- 
modyne detection setup. Secondly, for a fixed driving 
frequency, wl, the frequency content of the transmitted 
field can be determined by connecting the output of a 
photo-detector to a spectrum analyzer and searching for 
the broad-band, incoherent part as exemplified in Fig. [5] 
We remind that all calculations of the present paper 
are only valid in the linear regime, i.e. when atomic ex- 
citation is small, which can be confirmed theoretically 
by using Eq. (1201) . or experimentally by observing that 
transmission or reflection properties are independent of 
the intensity of the driving field, /3. 



VII. CONCLUSION 

For an optical cavity coupled to N two-level atoms, 
the steady-state reflection and transmission coefficients 
for both fields and intensities have been calculated ana- 
lytically as a function of the frequency, wl, of a coherent 
driving field, assuming small atomic saturation. In addi- 
tion, the frequency spectrum of the field emitted from the 
cavity has been determined. It has been demonstrated 
that atomic dephasing noise, independently subjected to 
each atom or common to all, prevents the cavity field 
from being in a coherent state, with two practical impli- 
cations: (1) intensities do not coincide with the modulus 
square of the field amplitudes, and (2) a broad-band, in- 
coherent part emerges in the cavity-emission spectrum. 
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Appendix A: Transforming a rapidly fluctuating 
Hamiltonian to a Lindblad phase damping term 

This appendix derives the Lindblad form of a phase 
decay, which is caused by a fast fluctuating term in the 
Hamiltonian. Consider the simple system: H = %e(t)0, 
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where e{t) is a real function fulfilling: (s(t)) — and 
(e(t)e(t')) = D5(t — t'), and O is a time-independent her- 
mitian operator. The time evolution operator is given by: 
U (ti, to) = exp(— id f t * e(t)dt), which in turn evolves the 

density matrix as: p(t\) = U{t\,to)p(to)W (t\,to). By 
expanding the time-evolution operator to second order 



(which makes sense when t\ — to is sufficiently small): 

pti i pt\ pti 

U(ti,t )f*il-id / e{t)dt~-6 2 I I e(t)e(t')dtdt' , 

(Al) 

it is possible to express the density matrix as: 



p(h)& p(t )-i[O,p(t )] I e(t)dt-- 



to 



ti ft! 



e(t)e{t')dtdt' 6 2 p(t ) + p{t )6 2 - 2Op(t )d 



to J to 



-> P(*o) 



D 



(A2) 



(ti - t ) O 2 p(*o) + P(^o)0 2 - 2Op(t )O 



where the arrow denotes averaging over the fast fluctua- 
tions of e(t). Now, by taking the limit, t\ — > to, we reach 
the standard Lindblad form of Eq. ([3]) with C = \f~DO. 



Appendix B: Derivation of (aja c ) and p c 



where the constants a, b, c, d, e, and / are given by: 



The steady-state solution of Eqs. (fTo^ - ljTg)) is derived in 
the following. Firstly, in Eq. (fTS)) set k — j and perform 
the summation over j. By adding the resulting right- 
hand side (which is zero in steady state) to Eq. (|16|) . we 
obtain: 

2k (a\a c ) + A 7 ||Pcxc = [P (fit) + /3* (a c )] , (Bl) 



a\dj 



is the mean atomic excita- 



where p cxc = i £\ 
tion. The above equation states the conservation of en- 
ergy: On the left-hand side the two terms describe the 
energy loss by cavity leakage and atomic population de- 
cay, which must be balanced by the right-hand side — the 
energy delivered by the coherent driving field. Secondly, 
an explicit expression for J2j k ^la,* ; can be obtained 
from Eq. (fT8|) : 



N 

£ 



4«j 



N Pc 



#7|| 



mi- 



(B2) 



This can be inserted into Eq. (fTTf leading to an expres- 
sion for J2j ("c%), which again is inserted into Eq. (fT8|) 
[being summed over j = ft]. As a result, another equa- 
tion involving (ajd c ), p cxc , and the already known cavity 
field (a c ), is obtained, and we are left with two equations 
involving two unknowns: 



a b 




Pcxc 




e 


c d 




_(aja c )_ 







<&c> I 2 , 



(B3) 



a = 7|| [(« + 7± )' + Af c ] + 2. 9 2 ( K + 7x )- 
-2 3 2 ( K + 7± ), 

2k, 

2 5 2 



1 + f 7||T 



71 + A2 



{ 5 ^( K + 7± ) + K ( 7 i-A 2 ) 



7± ( K 2 +A 2 )-2 7± A a A c }, 



7± + A 2 



(B4) 



The determinant, D = ad— be, of this set of equations was 
stated in Eq. (|2"Tj) . Now, (aja c ) and p cxc can be deduced 
using Kramer's rule leading to the results of Eqs. (fT9|) 
and (j2"0j) . 



Appendix C: Derivation of the cavity-field frequency 
spectrum 

We present here the computational details for the 
cavity-field frequency spectrum, Sa c (u>). Since it is 
a measure of the frequency distribution of the optical 
energy, an actual measurement of this entity can be 
modeled by coupling (weakly) an auxiliary narrow-band 
probe cavity to the optical cavity and computing the 
steady-state photon number of the probe cavity. In other 
words, the physical system described by the Hamiltonian, 
H, is extended to: 

H' = H + ftAa^ap + hg p (dld p + a c a v ), (CI) 

where A = to — lul denotes the observation frequency in 
the frame rotating at and g p measures the strength 
of the coupling. In addition, the (half-width-half- 
maximum) bandwidth, k p , of the probe cavity is mod- 
eled by an extra Lindblad term with C = 72K p a p . In 
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steady state we derive: 



-jg v (q c ) 

■Cn + iA 



(C2) 



which must be kept small for maintaining a negligible 
disturbance of the cavity field, d c . By defining g p = en p , 
the probe field is at least a factor of e weaker than d c 
for any value of k p and A and it suffices to calculate 
(d p d p ) to second order in e. Also, for a sufficiently small 
e, the probe field d p can be truncated to the two lowest 
Fock states {|0) p ,|l) p } with d p = |0) p (l| p , etc. The 

Hamiltonian H' and the entire density matrix p can then 
be expressed as: 



H' = 



'H 00 


£01- 






H 11 





and 



H 

hg p a c 



p oo p oi 
p w p n 



hg v a\ 
H + hAl 



(C3) 



(C4) 



ft follows that (d p d p ) = Tr(p n ) and (d p ) = Tr(p 10 ). 
Each of the four sub-blocks in H' have the dimensionality 
of the original, unperturbed system, H, and the evolution 
can be written: 



dp 



dt 

dp 10 
dt 



=-([tf°V° 



f 

ih 

£(p 00 

:i 

ih 

+ C(p 10 ) - k pP 



H°y° - p 01 H w ) 



,11 



--^(H w p 00 



2n p p 
- H n p 10 



(C5) 



i 10 H 00 



11/rlON 



(C6) 



where C( ■ ) is the Lindblad super-operator of the original 
system and the probe-cavity decay is added separately. 
Now, in steady state set Eq. (|C5|) equal to zero and com- 
pute the trace: 

<d p d p ) = T^p 11 ) = f [Tr(d c p 10 ) Tr(d c p 01 )] . (C7) 

The value of Tr(ajp 10 ) and its complex conjugate is re- 
quired and may be deduced using Eq. (|C6|) [in steady 



state, multiply from the left by dj and compute the trace] : 

^Tv([alH]p w ) + Tv(aiC(p 10 )) 
in 

= («p + iA)Tr(ay o ) + ig p (a\a c ) . (C8) 

This equation has a great similarity with the steady- 
state equation for However, we must replace 
(aj) — > Tr(djp 10 ), add an additional decay channel with 
rate k p , introduce another detuning of A, and the ex- 
tra term ig p (aja c ) must be included. The progress from 
here depends on the actual physical system as defined by 
H and £, but in practice the above-mentioned similari- 
ties allows for re-using existing equations in the spirit of 
the quantum regression theorem. 

Now, in order to keep the derivations simple we adopt 
the physical system of Sec. HH i.e. an empty, coherently 
driven cavity subjected to phase noise. In similarity with 
Eq. |gj), we obtain from Eq. (|C8l) : 



Tr(ajp iu ) 



-ig p (ata c ) + V2^P*Ty(p w ) 
k p + T — i(A c — A) 

<^c) + ^|(fic)| 2 " 



(C9) 



ftp + r — i(A c — A) 



where Eqs. ([6|) and (|C2p were used in the second step 
— note we only require Tr(ajp 10 ) to first order in e, and 
hence the unperturbed values of | (d c ) | 2 and (aja c y and 
also the relation ([7]) can be used. Now, insert the above 
result into Eq. (|C7[) and re-arrange the terms to obtain: 



<a P a p ) 



(ftp +r)/(7TKTjit) 

( Kp + r) 2 + (A c ^A) 2 



k p /tt 



A 2 



(CIO) 

The Lorentzian intensity-transmission profile of the 
probe cavity has an effective width [area-to-height ratio] 
of 7tk p . Furthermore, since e 2 connects the magnitudes of 
| (d c ) | 2 and | (d p ) | 2 through Eq. ([C2]). the left-hand side 
of the above expression can be interpreted as the amount 
of energy per frequency interval in the probe cavity when 
taking the limit k p — > 0, which leads to Eq. (|27p. The 
case of atoms [Eq. ([25]) ] can be derived in a similar fash- 
ion. 
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